Magneto-transport Signatures of a Magnetic Gap in Hybrid 
Topological Insulator-Ferromagnetic Insulator Heterostructure 

Devices 

A. Kandala, A. Richardella, D. W. Rench, D. M. Zhang^T. C. Flanagan, and N. Samartl^ 

Department of Physics and Materials Research Institute, 
The Pennsylvania State University, University Park, PA 16802, USA 

(Dated: December 7, 2012) 

Abstract 

We develop hybrid ferromagnetic insulator (GdN) /topological insulator (Bi2Se3) heterostructure 
devices to probe the effects of broken time reversal symmetry on electrical transport in the surface 
states of a three dimensional topological insulator. Low temperature (0.28 K < T < 30 K) mea- 
surements of electrical transport in laterally patterned devices show that interfacing Bi2Se3 with 
ferromagnetic GdN yields distinctive signatures in the quantum corrections to the longitudinal 
magnetoconductance, clearly differing from the standard weak anti-localization seen in control 
channels of Bi2Se3. The absence of electrical transport in the insulating ferromagnetic GdN over- 
layer allows us to attribute these signatures to contributions from surface state transport in Bi2Se3. 
We show that the observations are consistent with the opening of a magnetic gap in the surface 
state spectrum due to broken time-reversal symmetry. 
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Combining three dimensional (3D) topological insulators (TIs) with magnetism is of great 
current interest because of the unique effects predicted when the time reversal protected 
topological surface states are modified by symmetry breaking magnetic perturbations [T]. 
The resulting phenomena are as varied as the quantum anomalous Hall effect [2], the topo- 
logical magneto-electric effect [3] and magnetic monopoles[l] among others. Experiments 
aimed at observing such phenomena have primarily focused on magnetically doped TIs[5H9]. 
In such materials, broken time-reversal symmetry leads to changes in the symmetry class 
of quantum corrections to transport [SI El E], a hedgehog spin texture accompanied by a 
quenched Berry phase |B] and promising signatures of an emerging quantum anomalous Hall 
effect [7] . Complementary to the direct magnetic doping of TIs, several theoretical schemes 
have also proposed explorations of a different sample geometry wherein a patterned insu- 
lating ferromagnet (FM) is interfaced with a TI [T0] - [T5] . When the magnetic easy axis of 
the FM overlayer is out-of-plane, a gap opens in the surface states of the vicinal TI and 
this massive Dirac Hamiltonian should lead to chiral ID edge states along FM domain walls 
where the mass changes sign. If the easy axis of the FM lies in-plane, no gap is expected to 
first order but the moment can be levered out-of-plane by an external field, thereby opening 
or closing a gap. More recently, theoretical calculations [16j have demonstrated that hexag- 
onal warping effects can lead to a gap, even with an in-plane FM. Patterned FM overlayers 
could also be used to manipulate Klein tunneling in the Dirac surface states pO] and to 
create tunnel barriers in the surface state via proximity-induced exchange fields IH]. Such 
FM/TI heterostructure schemes have a clear advantage for electrical transport experiments: 
any unusual transport phenomena resulting from interfacing a TI with a FM can clearly 
be identified with surface states even if the chemical potential is large enough to allow 
for parallel bulk transport. In this Letter, we demonstrate an essential step towards such 
"magnetic gating" proposals by creating hybrid electrical transport devices wherein we in- 
terface an insulating FM (GdN) with a TI (Bi2Se3). By studying the quantum corrections to 
the magneto-resistance (MR) in patterned transport devices, we reveal the effect of broken 
time-reversal symmetry on the surface state and show that the data are consistent with the 
predicted opening of a magnetic gap. 

GdN is an insulating FM of interest in the context of low temperature spintronic devices, 
particularly since it can be deposited by reactive sputtering at ambient temperature [T71 HH]. 
This is of prime importance in studies of magnetic exchange coupling effects, as it minimizes 
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thermal diffusion of magnetic species. Changing the nitrogen composition allows access to 
metallic [ini [20], semiconducting|2I] and insulating [TH] regimes, opening possibilities for its 
use as a spin-injector and detector, a spin- filter [22] or a ferromagnetic gate in novel spintronic 
devices. The semiconducting and insulating states display a maximum ferromagnetic Curie 
temperature Tc~ 65 — 70 K[T8l EHl El], while more conducting N-deficient forms can have 
Tc well above 100 KfUlED]. 

Bi2Se3 thin films were grown by molecular beam epitaxy (MBE) on InP (lll)A sub- 
strates under typical growth conditions described elsewhere |23j. Following a brief exposure 
to ambient atmosphere, the samples were transferred to a Kurt Lesker CMS-18 system; 
after an in-situ Ar^ surface clean, GdN was deposited on the Bi2Se3 films by reactive rf 
sputtering of Gd in an Ar:N2 environment at ambient temperature. The GdN films were 
deposited at a rate of ~ 0.1 A/sec, with a sputtering power of 4.93 W/cm^, in a 15% N2:Ar 
gas environment at 5 mTorr pressure, with a source to substrate distance of 20 cm. The 
base pressure in the chamber prior to sputtering was ~ 4x10^^ Torr. The lack of ultra high 
vacuum is probably responsible for the lower Tc of our films. Since GdN oxidizes instantly 
on exposure to atmosphere, the films were capped in-situ with ~ 60 nm of Au deposited by 
dc sputtering. 

The structural and magnetic characterization of the heterostructures is shown in Fig. [1} 
We used cross-sectional high resolution transmission electron microscopy (HRTEM) to study 
the nature of the interface. The phase contrast image of Fig. [T]^a) reveals a polycrystalline 
GdN layer on top of a single crystal thin film of Bi2Se3. The selected area diffraction (SAD) 
pattern (Fig. [l|^b)) taken over the Bi2Se3 region shows the diffraction pattern of the InP 
substrate and additional spots consistent with Bi2Se3 single crystal grown epitaxially along 
the (111) direction of InP. A fast Fourier transform (FFT) from the GdN region confirms its 
polycrystalline nature (Fig. [l](c)). No additional layer that could indicate an interfacial re- 
action is discernible. X-ray energy dispersive spectroscopy (EDS) line profile measurements 
were carried out in scanning TEM (STEM) mode to study the elemental distribution along 
the growth direction of the film stack. Figure 1(d) is a bright field image in STEM mode 
with the red line and the arrow indicating the line profile position and scan direction. Thirty 
EDS spectra were collected with a step size of 1.9 nm and a nominal probe size of ~ 2 nm. 
The line profile is plotted in Fig. 1(e). Each data point was obtained by integrating the 
Ka and Kp lines without a standard for calibration; and hence the results are qualitative. 
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FIG. 1: (Color online) (a) HRTEM image of the heterostructure reveals a polycrystalline GdN 
thin film on a single crystal thin film of Bi2Se3 with no discernible interfacial layer, (b) SAD 
pattern depicts spots associated with the InP substrate and Bi2Se3 growing along the (111) direc- 
tion, (c) FFT from the GdN region confirms its polycrystallinity. (d) Bright field image of the 
heterostructure. The red line and the arrow indicate the extent of the EDS line scan and scan 
direction respectively, (e) EDS line profile plot depicts the elemental distribution along the growth 
direction. (f) M vsT measurements of the Bi2Se3/GdN heterostructure reveal Tc~ 13 K. (g) M vs 
fioH measurements at 1.8 K along the (112) and (111) directions display hysteresis, and confirm 
that the easy axis is in-plane. Inset: M vs fioH curves over a 2 T field range. 

Nevertheless, the elemental distribution is consistent with the film stack and shows relatively 
sharp interfaces. Magnetization measurements in a superconducting quantum interference 
device (SQUID) magnetometer then confirmed the ferromagnetism in the heterostructures. 
Measurements were taken along both in-plane and perpendicular-to-plane directions in a 5 
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FIG. 2: (Color online) (a) Optical microscope image of a typical device depicts a Hall bar of 
Bi2Se3 with a bare channel and a channel capped with GdN/Au, enabling comparison of mag- 
netotransport. (b) 2 probe I — V curve of the Bi2Se3/GdN/Au structure displays a strong non- 
linearity confirming the highly insulating nature of the GdN. (c) Temperature dependence of the 
2D conductivity (a) for the bare (red) and capped (blue) channels for T < 30 K is depicted in a 
semilog plot. Curves are offset for clarity. Dotted lines are fits to cr cx In T. (d) MR at T = 400 
mK over a 50 mT field range from capped (blue) and bare (red) channels. While the bare region 
shows a cusp-like positive MR typical of WAL, the region capped with ferromagnetic GdN has 
a superimposed negative MR. Inset: expanded view of MR over a 1 T field range reveals similar 
WAL at high fields for both capped and bare channels. 

mT measuring field after a 0.9 T field cool. The temperature dependence of the magneti- 
zation M reveals the Tq to be around 13 K, obtained by extrapolating the line of steepest 
slope of the curve (Fig. [l](d)). The hysteresis loops in the M v/s hqH plots (Fig. [l](e)) 
reveal the easy axis to lie in-plane. 

We measured magnetotransport in Hall bar devices of an 8 nm thick Bi2Se3 thin film, 
patterned into multiple channels of dimension 20 fim x 60 fim by standard photolithog- 
raphy and wet etching (K2Cr207:HBr:DI). The ferromagnetic GdN/Au gate electrode is 
then defined over one of the channels by another step of photolithography, sputtering and 
a lift-off process. The sputtering conditions were identical to those employed for the mag- 
netic characterization. Figure [2]^a) is an image of a typical device showing channels of bare 
Bi2Se3 and ones covered with GdN. This enables a comparison of the magneto-transport 
properties of the Bi2Se3 film in both regions to see the effects of the insulating ferromagnetic 
layer. Transport measurements were carried out in a ^He cryostat (operating temperature 
range 300 mK— 30 K) with a superconducting vector magnet with a maximum full sphere 
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field of 1 T. Two probe measurements tlirougli tlie Au-GdN-Bi2Se3 lieterostructure sliow 
a fiiglily non-linear I — V curve at 400 mK (Fig. ^h)) confirming the insulating nature 
of tlie GdN. Tlie asymmetry is indicative of two different Scfiottky barriers, one eacli at 
tlie Au/GdN and GdN/Bi2Se3 interfaces. This verifies that there are no contributions to 
transport from currents through the ferromagnet. The temperature dependence of the two 
dimensional (2D) sheet conductance (cr) of the bare and capped channels (Fig. |2](c)) shows 
a down-turn that follows a o" ~ ln(T) behavior consistent with electron-electron (e-e) inter- 
action in two dimensions |24j. Hall measurements reveal a carrier density of 1.36 x 10^^ cm~^ 
and 5.7 X 10^^ cm~'^ for the bare and capped channels. The larger carrier density observed 
in the capped channel may be due to interfacial defects. At 400 mK, both channels show a 
high-field positive MR attributed to weak anti-localization (WAL) due the large spin-orbit 
coupling in Bi2Se3 (inset. Fig. |2](d)). A closer look at the low-field regime reveals, however, 
a small negative MR below 10 mT in the capped channel, superimposed upon the WAL 
dip (Fig. [2]^d)). No such feature is seen for the bare Bi2Se3 channel. This behavior was 
reproduced in a second device with a similar carrier density for the capped channel. 

The temperature evolution of the change in MR of the capped channel is depicted in 
Fig. |3} The negative MR is very narrow at T = 400 mK but broadens continuously 
with increasing temperature and is seen to persist even at 30 K (the maximum operating 
temperature of the measurement system), well beyond the Tq of the GdN film. This indicates 
the importance of loose spins at the interface. The feature is always superimposed on top 
of a positive WAL and/or classical MR. These observations are emphasized in the line cuts 
in Fig. |3](b,c) of the raw MR at 5 K and 25 K. One possible interpretation would attribute 
the negative MR to spin disorder scattering from randomly oriented magnetic moments at 
zero-field, as observed in several dilute magnetic semiconductors [25]. The application of 
an external field aligns these moments, thereby reducing scattering and consequently, the 
resistance. In order to test this interpretation, we tilt the magnetic field towards the plane 
of the sample. The low-field MR is depicted for different tilt-angles in Fig. 3(d) and is 
seen to scale with the perpendicular component of the magnetic field (Fig. 3(e)). Such 
a scaling makes spin-scattering an unlikely origin for the negative MR and instead points 
towards weak localization phenomena (WL). Indeed, magnetically doped Bi2Se3 films O |9] 
and nanoribbons[6j have also displayed a similar negative MR that has been attributed to 
WL. 
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FIG. 3: (Color online) (a) MR {AR{H) = R{0)-R{H)) of the GdN capped channel as a function of 
temperature T and magnetic field fioH. Panels (b) and (c) show the MR line cuts corresponding to 
the dotted lines at T = 5 K and T = 25 K. (d) Angular dependence of the MR of the capped channel, 
as the field is tilted from perpendicular-to-plane (violet) to in-plane (red), in 10° steps offset for 
clarity, (e) The MR of the capped channel scales with the perpendicular-to-plane component of 
the magnetic field H±. 

We present support for attributing the low field MR to WL by fitting the different quan- 
tum corrections to a modified [26] Hikami-Larkin-Nagaoka[2^ (HLN) equation: 

In the above equation the a^'s and l^/s correspond to the pre-factors and coherence lengths 
respectively associated with WL and WAL, and Ib is the magnetic length. Although the 
fit employs four free parameters, fitting to the experimental data serves as a useful test of 
weak localization. Indeed, the fits work well in the temperature range below 15 K, as shown 
for the 4 K and 13 K data in Fig. |4](a). At higher T, increased relative contributions from 
the classical MC and/or the limited maximum field range precludes an effective fitting of 
the WAL. We also fit the regular HLN equation to the WAL data from the bare channel. A 
comparison of the different coherent lengths and pre-factors obtained is depicted in Fig. 4(b) 
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FIG. 4: (Color online) (a) Data (open circles) and modified HLN fits (bold lines) to the magneto- 
conductance of the capped channel at 4 K (blue) and 13 K (red). Curves are offset for clarity, (b) 
Temperature dependence of l,f, associated with WL in capped channel (red circles) , WAL in capped 
channel (blue circles) and WAL in bare channel (green circles). The dashed green line is a T~^'^^ fit 
to the high-temperature dependence of l^j, of the bare channel, suggestive of e-e interaction as the 
dominant dephasing mechanism, (c) Temperature dependence of a associated with WL in capped 
channel (red circles), WAL in the capped channel (blue circles) and WAL in bare channel (green 
circles). Error bars are 95% confidence intervals of the fits, (d) Magnetoconductance in Device B 
shows a weakened WAL from the capped channel (blue circles) in comparison to the bare channel 
(red circles), (e) T dependence of of the bare (red) and capped (blue) channels in Device B. 

and (c). The associated with the bare channel reaches a maximum of ~ 300 nm at 300 mK. 
The high temperature power law decay is suggestive of e-e interaction being the dominant 
dephasing mechanism. The associated with WL in the capped channel compare well with 
the bare channel, while those associated with WAL are seen to be significantly reduced. The 
a's associated with both channels do not show a sharp temperature dependence over the 
discussed temperature range. Consistent with existing ideas, the a associated with WL is 
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seen to be positive and fluctuates in the range ~ 0.1-0.2. However, the departure from the 
theoretical hmit (a = 1) may be attributed to signiflcant magnetic scattering driving the 
system towards the unitary class. The magnetoconductance of a second device B with a 
signiflcantly larger carrier density of ~ 11 x 10^^ cm^'^ for the GdN capped channel revealed 
no signature of WL. Figure 4(d) however depicts the signiflcantly weaker and broader WAL 
feature from the capped channel in comparison to the bare channel. Consequently, the 
associated with the capped channel in this device is also seen to be much shorter than the 
bare channel, possible due to increased defect scattering. 

Intuitively, one would expect magnetic scattering to destroy any form of quantum inter- 
ference. However, the opening of a sufficiently large gap in the surface state Dirac cone due 
to an out-of-plane magnetization at the surface is predicted ^26j to yield a crossover from 
WAL to WL. The crossover is expected to be dependent on the ratio {A/2Ef) of the gap 
size (A) to the Fermi energy with respect to the Dirac point (Ep). In the limit of small 
A/2Ep, the MR is positive and dominated by WAL, crossing over to WL in the limit of large 
A/2Ep via an intermediate regime of parabolic, positive MR. Such phenomena are predicted 
to arise due to a competition between magnetic scattering and the opening of a gap in the 
surface state gap. Our data can therefore be interpreted as evidence for the opening of such 
a surface state gap due to broken time-reversal symmetry at the interface. The absence 
of a crossover to WL in device B may then be attributed to its larger carrier density and 
hence a larger Ep. The higher chemical potential also implies a larger bulk conductivity 
that overwhelms any surface related signal. The WAL component that is always seen in 
these devices may arise from the bulk or from the bottom surface state, both of which are 
unlikely to be affected by the magnetic atoms at the top interface. 

Finally, we discuss the possible origins of the opening of a surface state gap in our 
devices. One possibility is that the application of a perpendicular magnetic field aligns the 
magnetization of the GdN out of plane. However, it is important to note that the WL is seen 
to persist beyond the Tq of the GdN, implying that the Dirac state gap also persists beyond 
Tc- Further, the values of a for the WL appear to be relatively temperature independent, 
implying that the size of the gap does not change much over the range of temperature studied. 
This raises the intriguing possibility of the development of an interfacial ferromagnetism 
with a much higher Tq than that of the GdN, that may be undetectable by conventional 
SQUID measurements. Indeed, recent mean field calculations predict a robust surface state 
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mediated ferromagnetism can arise due to RKKY interactions between magnetic atoms at 
the surface of a TI [2H] ■ Given the large E-p in these samples, the origin of such an interfacial 
ferromagnetism is unclear but is consistent with signatures of a high Tq magnetically ordered 
state revealed in recent measurements of x-ray magnetic circular dichroism in Mn-Bi2Se3 thin 
films [HI [9]. 

In summary, we have developed a new class of hybrid TI heterostructure devices wherein 
the surface transport in the TI is influenced by an overgrown insulating ferro magnet. By 
designing Bi2Se3 devices that contain channels selectively capped or not capped with GdN, 
we carried out a direct comparison of magnetotransport properties with and without mag- 
netic interactions. The crossover to WL in the capped channel is consistent with theoretical 
predictions that associate this behavior with the opening of a surface state gap. The devel- 
opment of these heterostructure devices opens up a broad range of opportunities for cleanly 
studying the influence of magnetism on electrical transport in the topological surface state. 
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